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Abstract 

In  this  Final  Report  we  present  results  for  both  areas  of  research  and  potential  applications  for 
the  novel  Left-Handed  MetaMaterials  (LHM)  which  our  group  has  previously  discovered.  In 
Section  I,  we  report  advances  made  in  our  laboratory  in  instrumentation  and  simulations  that 
focussed  on  the  utilization  of  LHM  that  can  serve  as  backplanes  for  antennas  with  improved 
radiative  properties.  In  Section  n  we  discuss  preliminary  measurements  on  the  incorporation  of 
regions  of  paramagnetic  spins  into  the  LHM  structure  so  as  to  permit  magnetic  field  modulation 
of  the  Split  Ring  Resonators  which  provide  the  negative  permeativity  of  the  sample.  In  Section 
m  we  report  on  our  efforts  to  create  a  cooperative  program  with  a  suggested  industrial 
manufacturer  for  utilization  of  high  Tc  superconducting  components  in  LHM  of  interest. 

I.  Metamaterials  as  Antenna  Substrates 

Our  goal  in  this  part  of  our  program  was  to  evaluate  metamaterials  for  application  as  novel 
antenna  substrates  or  radomes  for  communications  applications.  Ultimately,  our  goal  is  to 
incorporating  active  elements  into  a  metamaterial,  thereby  implementing  dynamic  beam  steering 
or  other  mechanisms  of  novel  RF  control.  Before  this  ultimate  goal  can  be  realized,  however, 
there  are  several  crucial  steps  that  need  to  be  completed  as  part  of  this  overall  program;  in 
particular, 

1.  The  ability  to  rapidly  simulate  and  fabricate  prototype  metamaterial  samples  must  be  in 
place. 

2.  We  need  to  have  an  efficient  fabrication  and  reliable  test  apparatus  to  enable  quick 
confirmation  of  the  expected  material  properties. 

3.  For  passive  materials,  we  need  to  evaluate  metamaterial  structures  that  indicate  particular 
advantages  and  capabilities  either  not  present  in  existing  materials  or  exhibit  substantial 
improvements. 

4.  For  dynamic  tuning,  we  need  to  evaluate  two  or  more  material  properties — or  “states”  of 
material,  which  will  significantly  modify  the  radiation  pattern  associated  with  a  nearby 
antenna,  and  the  means  to  switch  between  these  states. 

5.  For  the  most  flexible  control  uses,  we  need  to  perform  a  systematic  analysis  of  the 
possibilities  inherent  to  a  reconfigurable  material  whose  material  properties  can  be  adjusted 
spatially,  on  a  scale  less  then  the  wavelength  (i.e.,  at  the  unit  cell  level). 


In  this  phase  of  the  project,  we  have  developed  the  groundwork  that  will  enable  the  future  goals 
of  the  project  to  be  realized.  Specifically,  in  this  project  period  we  have: 

1.  Improved  our  metamaterial  development  cycle  time.  Under  separate  fiinding,  we  have 
obtained  an  LPKF  ProtoMat  circuit  board  prototype  numerically  controlled  milling  machine, 
capable  of  producing  appropriate  sized  metamaterial  samples  within  one  or  two  days. 

2.  Modified  our  Angle-Resolved  Microwave  Spectrometer  (ARMS)  to  allow  antenna  gain 
measurements  from  a  coaxial  fed  antenna.  We  have  systematically  improved  the 
reproducibility  of  the  gain  measurements  to  a  level  of  less  than  1%  from  scan  to  scan;  with 
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this  sensitivity  we  csn  easily  evaluate  the  antenna  enhancements  and  even  slight  changes  in 
gain. 

3.  Performed  calibration  measurements  on  low-impedance  reflector  planes  to  compare  with 
expected  gain  curves. 

4.  Measured  metamaterial  reflector  planes  and  compared  with  expected  gain  curves  based  on 
the  calculated  material  properties. 


The  measurements  are  preliminary  and  primarily  made  to  test  the  feasibility  of  our  methodolo^. 
In  this  section  of  the  report,  we  describe  the  modifications  to  the  ARMS  apparatus,  its  sensitivity 
and  limitations  and  the  preliminary  metamaterial  measurements. 


Modifications  to  the  ARMS  Apparatus 

The  ARMS  apparatus  was  originally  developed  to  better  quantify  the  refraction  properties  of 
left-handed  metamaterials.  However,  the  apparatus  is  capable  of  probing  a  number  of  different 
metamaterial  aspects,  including  the  fiinction  of  a  metamaterial  used  as  an  antenna  back  plane. 

The  ARMS  design  follows  in  spirit  that  utilized  originally  to  confirm  the  negative 
refractive  index  of  a  wedge-shaped  metamaterial  sample.  The  ARMS  apparatus  is  based  on  the 
planar  waveguide  geometry  shown  in  Fig.  1.  Two  parallel  plates  spaced  nominally  one 
centimeter  apart  are  used  to  confine  radiation  to  the  plane.  By  restricting  radiation  in  this 
manner,  we  minimize  the  necessary  sample  volume  required  for  the  measurements.  The  lowest 
mode  of  the  planar  waveguide  has  no  cutoff  and  is  characterized  by  an  electric  field  polarized 
uniformly  along  the  direction  perpendicular  to  the  plane  of  the  waveguide  plates.  In  the  empty 
chamber,  these  modes  have  the  identical  dispersion  characteristics  of  free  space  waves. 
Moreover,  because  the  electric  field  is  always  polarized  perpendicular  to  the  plane  of  incidence, 
as  defined  by  the  propagation  direction  and  the  normal  to  any  sample  surface,  the  waveguide 
modes  are  equivalent  to  s-polarized  waves.  Note  that  while  the  confined  geometry  allows 
minimum  sample  sizes  to  be  considered,  the  restriction  to  one  polarization — s-polarizatiqn  ^is  a 
distinct  limitation  that  prevents  a  fiill  sample  characterization.  However,  the  efficiency 
associated  with  smaller  sample  sizes  and  rapid  sample  characterization  offset  this  disadvantage. 

Of  particular  concern  for  this  work  is  the  possibility  of  reflection  of  the  generated  and 
scattered  waves  back  into  the  scattering  region,  which  could  produce  unwanted  artifacts  in  the 
data.  At  all  interfaces  between  the  scattering  chamber  and  free  space,  for  example,  there  exists  a 
substantial  impedance  mismatch  that  will  reflect  some  portion  of  the  radiation.  A  possible 
means  of  minimizing  or  eliminating  this  reflection  would  be  to  line  the  entire  periphery  of  the 
chamber  with  microwave  absorber;  but  this  would  reduce  the  signal  reaching  the  external 
detector,  complicate  the  design  and  possibly  also  produce  artifacts  in  the  data.  A  compromise 
absorber  pattern  is  utilized  in  the  present  apparatus,  as  indicated  in  Fig.  1  (right  panel),  in  which 
a  toothed  absorber  structure  is  arranged  on  the  back  plane  of  the  central  chamber.  As  the 
subsequent  results  show,  the  absence  of  absorber  along  the  arc  of  the  semi-circular  plates  does 
not  appear  to  interfere  noticeably  with  the  data,  e.g.,  standing  wave  modes  within  the  chamber 
are  not  problematic  as  the  large  chamber  size  and  the  absorber  along  the  backside  sufficiently 
reduce  these  unwanted  artifacts  compared  to  the  expected  radiation  patterns. 

The  radius  of  the  circular  plates  sets  the  detector  at  16”  (chamber  radius)  from  the 
radiating  element.  The  combination  of  the  coaxial  antenna  and  the  substrate  set  an  effective 
aperture,  which  in  turn  determines  the  characteristic  Fresnel  length  separating  the  propagating 
near-field  region  from  the  far-field  region  of  interest.  We  have  not  performed  detailed 


2 


calculations  to  determine  the  effective  aperture  of  the  antenna,  but  rough  arguments  suggest  that 
we  are  safely  in  the  far-held  region  where  the  fields  are  measured.  The  data  thus  far  acquired 
confirms  this,  but  we  will  consider  the  topic  in  more  detail  in  future  work. 

The  waveguide  detector  is  mounted  flush  with  the  chamber  (on  the  semicircular  face)  on 
a  toothed  gear  stepper  motor,  driven  by  customized  software  that  also  controls  the  data 
acquisition.  Microwaves  are  injected  into  the  chamber  via  a  removable  coax-fed,  source  antenna 
located  at  the  center  of  the  circular  (Fig.  1,  right  panel).  When  using  the  co^ial  antenna,  we 
extend  the  back  surface  absorber  material  across  the  plane  wave  aperture  and  insert  the  antenna 
through  a  hole  in  the  bottom  plate. 


Figure  1:  (Left  panel)  A  schematic  of  the  ARMS  apparatus,  in  which  a  plane  wave,  launched 
down  the  rectangular  wave-guide,  interacts  with  a  meta-material  sample  placed  at  the  center  of 
the  circular  plates  (A=16”,  B=6”,  C=16”,  D=8”,  plate  spacing=  0.4”).  The  detector,  confined  to 
the  radius  by  a  bearing  mounted  arm,  sweeps  around  circumference  of  the  plates  driven  by  a 
stepper  motor.  Toothed  absorber  (in  black)  is  placed  inside  the  plates  to  minimize  unwanted 
reflective  scattering.  (Right  panel)  ARMS  apparatus  using  the  optional  antenna  source  (A=1.08”, 
B«0.315”,  C=0.375”,  D=  0.162”  dia,  E=  0.050”  dia),  combined  with  experimental  backplane 
sample.  The  angle-resolved  radiation  pattern  is  determined  by  scanning  the  circumference  with 
the  detector. 

Microwaves  are  generated  and  detected  in  the  ARMS  apparatus  by  an  Agilent  Vector 
Network  Analyzer  (VNA,  model  number  8722).  While  the  frequency  range  of  the  VNA  is  50 
MHz  through  40  GHz,  the  plate  separation  of  the  planar  waveguide  restricts  our  measurements  to 
the  frequency  band  8-12  GHz,  although  we  find  acceptable  data  can  be  acquired  in  the  range  8- 
16  GHz.  A  data  acquisition  program,  written  as  a  series  of  macros  in  Excel  VBA  (Visual  Basic 
for  Applications),  controls  the  stepper  motor  (Little  Step-U,  TLA  Microsystems,  Auckland,  New 
Zealand)  that  rotates  the  detector  arm,  in  addition  to  controlling  the  VNA  (Fig.  2). 
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The  angular  stepping  motion  of  the  detector  induces  vibrations  in  the  detector  arm.  In 
order  to  avoid  introduction  of  unwanted  structure  into  the  data,  a  delay  (usually  one  second)  is 
used  to  allow  for  damping  of  the  vibrations  before  the  control  program  initiates  acquisition  or 
averaging  of  the  data. 

The  accuracy  of  the  ARMS  waveguide  apparatus  is  in  part  determined  by  the  flatness  of 
the  Vi”  thick  aluminum  plates  (rolled  sheet  aluminum)  which  make  up  the  upper  and  lower 
central  chamber  plates  and  the  plates  of  the  extended  channel.  As  purchased,  the  lower  central 
chamber  plate  was  significantly  warped.  To  remedy  this,  it  was  necessary  to  flatten  the  plate  by 
mounting  three  1”  square  aluminum  straightening  rods  to  the  bottom  surface.  The  rods  were 
drilled  with  several  >  Vi”  through-holes,  with  corresponding  Vi  x  20  blind  tapped  holes  machined 
into  the  bottom  of  the  lower  aluminum  plate.  When  the  straightening  rods  were  brought  flush 
with  the  aluminum  plate  surface  by  means  of  tightening  the  Vi  x  20  machine  screws,  the 
aluminum  plate  became  nearly  flat.  Starting  with  a  machine  ground  Vi”  aluminum  plate  rather 
than  the  rolled  sheet  aluminum  might  reduce  the  need  for  straightening  bars  for  further  improved 
performance. 

The  antenna  source,  which  in  our  case  is  the  projecting  pin  of  a  coaxial  panel  jack 
receptacle  (Applied  Engineering  Products,  Model  9004-xl  13-000  or  similar)  is  used  to  produce 
cylindrical  waves.  In  keeping  with  the  two-dimensional  confined  nature  of  the  ARMS  apparatus, 
we  would  ideally  like  to  use  a  linear  current  source  as  the  localized  excitation,  which  would 
couple  strictly  to  the  lowest  propagating  modes  in  the  chamber.  In  reality,  our  source  protrudes 
into  the  chamber  a  finite  distance,  and  produces  image  charges  on  the  upper  and  lower  plates. 
The  resulting  electromagnetic  source  presumably  couples  in  a  complicated  manner  to  the  higher 
order  modes  of  the  chamber,  although  none  of  these  higher  order  modes  propagate  at  frequencies 
below  14.8  GHz.  Nevertheless,  the  antenna  excitation  we  use  has  excellent  reproducibility,  and 
can  be  used  to  probe  antenna/metamaterial  interactions. 

Fig.  2  shows  the  reproducibility  of  an  angular  scan,  with  the  power  introduced  by  the 
antenna  excitation.  Two  scans  were  made,  and  the  intensity  from  one  scan  divided  by  the 
intensity  of  the  second  scan.  The  reproducibility,  as  determined  from  the  scatter  in  the 
measurement,  corresponds  to  —30  dB.  For  antenna  backplane  measurements,  where  ~0. 1  dB 
accuracy  is  typically  needed,  the  ARMS  chamber  can  provide  extremely  sensitive  measurements. 


Antenna  Backplcme  Calculations 

The  solution  of  a  line  source  above  an  infinite  half-space  of  arbitrary  material  can  be  readily 
found  by  expanding  the  line  source  in  terms  of  plane  waves,  and  calculating  the  reflection 
coefficient  for  each  plane  wave.  We  have  developed  this  semi-anal54:ical  method  [see  W.  C. 
Chew,  Waves  and  Fields  in  Inhomogeneous  Media,  (IEEE  Press,  New  York,  1995),  Section 
2.2.1]  in  order  to  compare  with  experiments  and  numerical  finite-element  based  simulations. 

An  example  of  the  utility  of  the  method  is  shown  in  Fig.  3a,  where  we  plot  the  intensity 
versus  angle  of  a  line  source  next  to  a  “perfect  electric  conductor,”  or  PEC.  The  PEC  tends  to 
short  the  line  source  as  it  approaches  the  surface,  so  that  the  total  radiated  power  is  reduced. 
When  the  line  source  is  ^74  away  from  the  surface,  however,  the  reflected  waves  are  constructive 
so  that  the  forward  gain  is  actually  enhanced  relative  to  the  line  source  in  free  space. 
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Figure  2:  Experimental  data  which  displays  the  reproducible,  low  noise  angular  scan  of  the 
ARMS  apparatus  with  the  center  antenna  source  radiating  without  a  sample  backplane.  The  low 
noise  ratio  of  two  consecutive  scans  means  that  by  normalizing  the  instrument  to  an  empty  scan, 
the  radiation  pattern  from  a  sample  backplane  can  be  accurately  measured. 

The  influence  of  material  properties  can  be  seen  by  comparing  Fig.  3a  with  Fig.  3b, 
which  shows  gain  curves  for  a  line  source  next  to  a  “perfect  magnetic  conductor,”  or  PMC.  The 
PMC  is  not  realizable  in  conventional  materials,  but  provides  an  indication  as  to  how  materials 
might  be  engineered  for  useful  communications  applications.  Note  that  instead  of  reducing  the 
radiative  power  when  the  line  source  approaches  the  surface,  the  PMC  actually  enhances  the 
radiative  power  by  roughly  a  factor  of  two  over  a  free  space  line  source;  furthermore  the  gain 
curve  is  flat  over  all  angles  when  the  line  source  is  in  direct  contact  with  the  surface.  This 
example  is  illustrative,  as  it  is  not  our  intent  to  pursue  PMC  equivalent  materials;  numerous 
groups  are  working  on  this  application,  with  some  success.  Our  immediate  purpose  is  to  utilize 
the  method  to  compare  with  measurements  in  the  ARMS  apparatus.  We  present  these 
measurements  in  the  following  section.  Because  we  will  eventually  be  interested  in  the  response 
from  different  types  of  metamaterial  substrates,  we  have  created  a  general  code  (in  Mathematica) 
to  calculate  gain  curves  for  a  line  source  above  a  half-space  of  material  characterized  by  arbitrary 
E  and  n  tensors.  Thus,  we  have  the  capability  to  study  materials  such  as  the  indefinite  media 
discussed  in  previous  reports. 
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Figure  3:  (a)  Gain  curves  for  a  line  source  at  varying  distances  from  a  “perfect  electric 
conductor”  (PEC),  (b)  Gain  curves  for  a  line  source  at  varying  distances  from  a  “perfect 
magnetic  conductor”  ^MC). 


Calibration  Measurements:  PEC 

To  determine  whether  our  ARMS  apparatus  would  produce  reliable  antenna  back  plane 
measurements,  we  first  studied  the  coaxial  line  at  varying  distances  (and  frequencies)  in  front  of 
an  approximate  PEC.  In  this  case,  the  PEC  was  a  rectangular  aluminum  slab  placed  inside  the 
chamber  at  a  fixed  distance  of  1.3  cm  from  the  antenna.  Rather  than  moving  the  antenna  to  study 
the  distance  dependence  of  the  gain  patterns,  the  frequency  was  swept  from  8  to  16  GHz.  The 
values  of  6  shown  in  the  figure  refer  to  the  distance  of  the  line  source  to  the  PEC  in  terms  of 
wavelength.  The  three  plots  in  Fig.  4  demonstrate  quite  good  agreement  with  the  semi-analytical 
model  described  above.  There  is  some  deviation  from  the  theoretical  curves,  especially  near  the 
end-fire  positions;  this  deviation  may  be  expected  as  our  sample  is  of  finite  extent  and  as  the 
antenna  we  are  using  is  significantly  different  from  an  ideal  line  source.  Given  the  simplicity  of 
the  model  and  the  actual  experimental  configuration,  the  agreement  is  quite  acceptable. 


Figure  4:  Gain  profiles  for  a  coaxial  antenna  backed  by  a  conducting  plane.  The  dots 
correspond  to  the  measured  gain,  while  the  overlayed  solid  curves  correspond  to  the  gain 
calculated  by  an  analytical  model.  (Left)  Frequency  of  11.0  GHz,  (Middle)  14.3  GHz,  and 
(Right)  9.9  GHz. 
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Preliminary  Measurements  on  a  Metamaterial  Sample  asBaclplane 

Having  shown  the  feasibility  of  backplane  measurements  using  our  ARMS  apparatus,  our  next 
step  is  to  characterize  metamaterial  samples  as  anteima  backplanes.  This  step  requires  a 
considerable  effort,  as  the  material  properties  of  metamaterials  themselves  are  difficult  to 
characterize,  which  makes  comparison  with  theory  an  even  greater  challenge.  Nevertheless,  as 
our  final  step  in  this  reporting  period,  we  have  initial  results  that  indeed  appear  very  promising. 

One  of  the  simpler  structures  for  us  to  consider  is  a  metamaterial  sample  made  only  of 
split  ring  resonators  (SRR).  Using  our  established  methods  of  analysis,  we  simulated  an 
metamaterial  composed  of  SRRs.  The  design  consisted  of  SRRs  in  one  dimension  only,  so  that 
the  sample  was  anisotropic.  The  effective  eand  p  tensors  of  the  effective  medium  were 
recovered  using  S-parameter  inversion,  as  well  as  the  effective  index-of-reffaction  (n)  and 
impedance  (z).  The  latter  two  parameters  provide  the  best  indicators  of  expected  performance, 
particularly  for  isotropic  samples.  For  isotropic  samples,  either  large  or  small  z  implies  a  highly 
reflective  sample  that  can  clearly  function  as  a  backplane.  When  z  is  of  moderate  value, 
however,  the  index  provides  a  measure  of  the  penetration  depth  and  an  indication  as  to  how  thick 
the  sample  must  be  to  perform  as  a  reflector  plane.  The  SRR  sample  used  here  was  not  designed 
for  any  particular  application  other  than  as  a  test  structure;  our  goal  was  to  understand  how 
closely  the  gain  pattern  from  an  antenna  next  to  a  fabricated  metamaterial  would  match  the  gain 
predicted  fi-om  the  derived  material  parameters  used  in  the  semi-analytical  model. 

The  results  of  the  measurements  are  shown  in  Fig.  5,  in  wWch  we  present  plots  of  the 
predicted  material  properties  (e  and  p),  complex  index  and  impedance  and  measured  gain  curves 
for  the  SRR  metamaterial.  We  have  selected  three  fi-equencies  at  which  data  is  presented.  As  a 
fimction  of  fi'equency,  the  material  parameters  trace  out  curves  such  as  that  shown  in  the  lefl:most 
two  plots  in  the  figure.  Two  black  dots  indicate  the  calculated  values  of  either  e  and  p  (far  left) 
or  n  and  z  (middle  plot);  using  these  values  in  the  semi-analytical  model  of  a  line  source  next  to  a 
half-space  we  can  derive  the  expected  gain  curve,  shown  as  the  solid  green  curve  in  Fig.  5.  The 
measured  points  are  also  shown  in  the  figure.  These  curves  provide  strong  evidence  that  we  can 
apply  metamaterial  as  backplanes  and  obtain  reliable  and  predictable  results.  We  should  caution 
that  we  have  presented  three  of  the  better  regions  of  agreement;  at  other  frequencies  we  noted  a 
greater  discrepancy  between  predicted  and  measured  gain. 

Having  established  our  capabilities  of  metamaterial  design,  fabrication  and  measurement 
for  the  purpose  of  antenna  backplane  applications,  we  plan  to  continue  to  on  the  path  outlined  at 
the  beginning  of  this  section.  In  particular,  we  will  work  to  improve  further  the  measured  gain 
characteristics  with  that  expected  from  the  recovered  material  parameters  from  the  metamaterial 
designs.  We  will  evaluate  the  effects  of  finite-sized  backplanes  utilizing  full-wave  finite-element 
modeling.  We  will  also  analyze  the  gain  characteristics  for  backplanes  composed  of  materials 
whose  material  parameters  vary  throughout  the  volume.  This  effort  is  an  important  precursor  to 
advanced  reconjigurable  metamaterials. 
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Figure  5:  (left)  the  value  of  |i  as  a  function  of  frequency.  The  vertical  axis  corresponds  to  Im(p) 
while  the  horizontal  axis  corresponds  to  Re(p);  s  is  approximated  to  have  the  value  2.5  +  Oi 
throughout.  As  a  function  of  frequency,  |x  traces  out  the  curve  shown  in  the  complex  plane;  the 
two  black  dots  indicate  the  values  of  e  and  p,  at  the  frequency  studied,  (center)  We  present  the 
values  of  complex  n  and  z  as  a  function  of  frequency,  (right)  Measured  (dots)  and  calculated 
(solid  curve)  gain  curves  for  three  frequencies  (13.14  GHz,  13.82  GEtz  and  11.34  GHz).  The 
lighter  dots  of  the  measured  points  indicate  where  the  penetration  depth  is  larger  than  our  sample 
size,  such  that  these  points  are  likely  to  have  significant  error. 
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II.  Magnetic  Field  Modulation  of  the  Split  Ring  Resonators 

It  is  clear  that  as  negative  index  of  refraction  or  indefinite  media  designs  become  available  for 
applications  as  lenses,  antenna  substrates,  “designer”  index  and  surface  impedance  materials  and 
other  purposes,  there  will  be  a  corresponding  need  to  provide  the  features  of  signal  modulation, 
on/ofF  switching,  directional  radiation  control,  etc.  We  have  suggested  that  since  the  key 
properties  of  LHM  depend  on  the  ability  to  incorporate  specific  values  (of  both  positive  and 
negative  sign)  of  the  effective  permeability  tensor,  it  is  an  optimum  strategy  to  use  magnetic 
fields  to  provide  the  external  control.  One  of  the  simplest  approaches  is  to  place  strongly 
paramagnetic  spin  materials  at  the  centers  of  the  SRR  where  their  microwave  fields  are  the 
strongest.  While  all  types  of  magnetic  materials  may  eventually  prove  useful,  we  have  decided  to 
initially  use  paramagnetic  spins  as  being  the  most  instructive  and  flexible  first  case  to 
investigate. 

We  also  have  chosen  to  exploit  electron  spin  resonance  for  the  simple  reason  that  this  happens  to 
be  an  area  in  which  we  have  a  complete  range  of  spectrometers,  ancillary  equipment,  and  many 
years  of  direct  experience  designing  diverse  test  structures  and  measurement  techniques.  In  Fig. 
6  we  present  four  pictures  illustrating  our  initial  work.  The  large  left  side  picture  is  a  snapshot  of 
one  of  the  several  magnet-electronic-cryogenic  ESR.  stations.  The  top  right  side  figure  shows  the 
end  of  a  pair  of  waveguides  which  permit  measurements  of  the  transmission  of  microwaves  (in 
the  8-12  GHz  band)  through  a  small  sample  placed  inside.  In  the  bottom  smaller  picture  we  see  a 
blow  up  of  a  single  SRR  which  is  built  as  an  inner  and  outer  pmr  of  “C”  sections.  In  the  right 
small  picture  we  see  a  similar  SRR  coated  with  an  appropriate  paramagnetic  spin  material. 


Figure  6:  (left)  Magnet  and  microwave  spectrometer  apparatus,  (upper  right)  X-band 
microwave  waveguide  fittings  for  transmission  measurements,  (lower  middle)  A  split  ring 
resonator,  (lower  right)  A  split  ring  resonator  loaded  with  DPPH. 

In  Fig.  7  we  observe  the  transmitted  power  (curved  white  line)  as  a  function  of  microwave 
frequency  which  is  transmitted  past  an  individual  SRR  appropriately  oriented  within  the  center 
region  of  the  waveguide.  The  white  dot  at  the  center  of  the  lowest  point  of  curvature  of  the 
transmitted  power  indicates  the  center  frequency  of  the  actual  resonance  response  of  the  SRR.  In 
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We  study  the  behavior  of  wave  propagation  in  materials  for  which  not  all  of  the  principal  elements  of 
the  permeability  and  permittivity  tensors  have  the  same  sign.  We  find  that  a  wide  variety  of  effects  can 
be  realized  in  such  media,  including  negative  refraction,  near-field  focusing,  and  high  impedance 
surface  reflection.  In  particular,  a  bi layer  of  these  materials  can  transfer  a  field  distribution  from  one 
side  to  the  other,  including  near  fields,  without  requiring  internal  exponentially  growing  waves. 
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The  range  of  available  electromagnetic  material 
properties  has  been  broadened  by  recent  developments 
in  structured  media,  notably  photonic  band  gap  materials 
and  metamaterials.  These  media  have  allowed  the 
realization  of  solutions  to  Maxwell’s  equations  not  avail¬ 
able  in  naturally  occurring  materials,  fueling  the  discov¬ 
ery  of  new  physical  phenomena  and  the  development  of 
devices. 

Photonic  crystal  effects  typically  occur  when  the 
wavelength  is  on  the  same  order  or  smaller  than  the  lattice 
constant  of  the  crystal.  Metamaterials,  on  the  other  hand, 
have  unit  cell  dimensions  much  smaller  than  the  wave¬ 
length  of  interest  A  homogenization  process,  not  unlike 
that  applied  in  Ref.  [1],  can  allow  the  otherwise  compli¬ 
cated  composite  medium  to  be  described  conveniently  by 
a  permittivity  tensor  (e)  and  a  permeability  tensor  (pt), 
rather  than  by  band  diagrams. 

In  2000,  it  was  shown  experimentally  that  a  metama¬ 
terial  composed  of  periodically  positioned  scattering  ele¬ 
ments,  all  conductors,  could  be  interpreted  as  having 
simultaneously  a  negative  effective  e  and  a  negative  ef¬ 
fective  fjL  [2].  A  medium  with  simultaneously  isotropic 
and  negative  e  and  jjl  supports  propagating  solutions 
whose  phase  and  group  velocities  are  antiparallel;  equiva¬ 
lently,  such  a  material  can  be  rigorously  described  as 
having  a  negative  index  of  refraction  [3,4].  An  experi¬ 
mental  observation  of  negative  refraction  was  reported 
using  a  metamaterial  composed  of  wires  and  split  ring 
resonators  deposited  lithographically  on  circuit  board 
material  [5]. 

The  prospect  of  negative  refractive  materials  has  gen¬ 
erated  considerable  interest,  as  this  simply  stated  material 
condition  suggests  the  possibility  of  extraordinary  wave 
propagation  phenomena,  including  near-field  focusing  [6]. 
So  remarkable  have  been  the  claims  surrounding  negative 
refraction,  that  some  researchers  have  been  prompted  to 
examine  critically  the  achievability  of  negative  refraction 
in  existing  metamaterials  [7,8].  While  such  concerns 
might  appear  relevant  in  the  context  of  frequency- 
dispersive  materials,  the  interpretation  of  these  structured 
materials  as  negative  refractive  has  been  entirely  consis- 


PACS  numbers:  78.20.Ci,  42.25.-p,  78.67.-n,  84.40.-x 

tent  with  all  aspects  of  reported  experimental  data  [2,5], 
as  well  as  with  numerical  simulations  of  both  mono¬ 
chromatic  and  modulated  beams  [9,10].  For  the  purposes 
of  this  Letter,  we  thus  assume  that  the  descriptions  pre¬ 
sented  here  will  be  applicable  to  real  materials,  although 
the  extent  to  which  this  is  true  remains  a  topic  of  active 
pursuit  [11]. 

Lindell  et  aL  [12]  have  shown  that  the  property  of 
negative  refraction  is  not  confined  to  materials  with  nega¬ 
tive  definite  e  and  /x  tensors,  but  can  be  expected  to  occur 
in  certain  classes  of  uniaxially  anisotropic  media 
We  explore  here  the  general  electromagnetic  solutions 
associated  with  these  anisotropic  media,  focusing  on 
combinations  of  such  media  that  provide  access  to  the 
underlying  properties.  While  there  is  a  substantial  reflec¬ 
tion  coefficient  for  non-normal  waves  incident  on  an 
anisotropic  material  from  an  isotropic  material  (or  free 
space),  due  to  the  difference  in  dispersion  properties  of 
the  two  materials,  a  highly  transmissive  composite  struc¬ 
ture  composed  of  layers  of  positively  and  negatively 
refracting  anisotropic  materials  can  be  formed.  As  a 
specific  example,  we  present  an  analysis  of  a  compen¬ 
sated  bilayer  that  produces  subwavelength  near-field  fo¬ 
cusing,  but  mitigates  the  thickness  and  loss  limitations  of 
the  isotropic  “perfect  lens.” 

To  simplify  the  proceeding  analysis,  we  assume  a 
linear  material  with  e  and  pt  tensors  that  are  simulta¬ 
neously  diagonalizable, 

/£^  0  0\  //X,  0  0\ 

8=0  0  ,  =  0  0  .  (1) 

Vo  0  £,/  Vo  0  /X,/ 

Metamaterials  can  be  readily  constructed  that  closely 
approximate  these  e  and  |x  tensors,  with  elements  of 
either  algebraic  sign.  In  fact,  the  scattering  elements 
comprising  the  metamaterials  used  to  demonstrate  nega¬ 
tive  refraction  [5,13]  are  appropriate  building  blocks  for 
the  classes  of  materials  to  be  discussed  here  (see  Fig.  1). 

We  apply  the  term  indefinite  to  anisotropic  media  in 
which  not  all  of  the  principal  components  of  the  e  and  |x 


077405-1 


0031-9007/03 /90(7)/077405(4)$20.00  ©  2003  The  American  Physical  Society 


077405-1 


Volume  90,  Number  7 


PHYSICAL  REVIEW  LETTERS 


week  ending 
21  FEBRUARY  2003 


FIG.  1.  A  bilayer  composed  of  indefinite  media  implemented 
using  split  ring  resonators  and  straight  wires.  The  structures  in 
the  top  of  the  figure  implement  never  cutoff  mcdm  for  electric  y 
polarization.  Top  left  is  negative  refracting  and  top  right  is 
positive  refracting.  The  structures  in  the  bottom  of  the  figure  do 
the  same  for  magnetic  y  polarization. 


tensors  have  the  same  sign.  We  will  consider  here  layered 
media  with  surfaces  normal  to  one  of  the  principal  axes, 
which  we  define  to  be  the  z  axis.  We  assume  in  our 
analysis  a  plane  wave  with  the  electric  field  polarized 
along  the  y  axis  having  the  specific  form 

E  =  (2) 

though  it  is  generally  possible  to  construct  media  that  are 
polarization  independent,  or  that  exhibit  different  classes 
of  behavior  for  different  polarizations.  The  plane  wave 
solutions  specified  by  Eq.  2  satisfy  the  dispersion  relation 

(3) 

Since  we  have  no  x  or  y  oriented  boundaries  or  interfaces, 
real  exponential  solutions,  which  result  in  field  divergence 
when  unbounded,  are  not  allowed  in  those  directions. 
Furthermore,  is  restricted  to  be  real.  Also,  since 
represents  a  variation  transverse  to  the  surfaces  of 
our  layered  media,  it  is  conserved  across  the  layers,  and 
naturally  parametrizes  the  solutions. 

In  the  absence  of  losses,  the  sign  of  kl  can  be  used  to 
distinguish  the  nature  of  the  plane  wave  solutions.  >  0 
corresponds  to  real  valued  and  propagating  solutions. 
kl<0  corresponds  to  imaginary  and  exponentially 
growing  or  decaying  (evanescent)  solutions.  When 
Syfiz  >  0,  there  will  be  a  value  of  for  which  kl  =  0, 
This  value,  which  we  denote  is  the  cutoff  wave  vector 
separating  propagating  from  evanescent  solutions.  From 


Eq.  (3),  this  value  is  kc  =  We  identify  four 

classes  of  media  based  on  their  cutoff  properties: 


Media  conditions  Propagation 


Cutoff  Sy/Jijf  >  0 

Anti-cutoff  Syfix<0 

Never  cutoff  ^yf^x  ^  ^ 

Always  cutoff  <  0 


>0  k^<kc 

fix/fiz  <0  kj,>  kc 

Mjc/ ^  0  All  real 

J^x/f^z  >0  No  real 


Note  the  analysis  presented  here  is  carried  out  at  con¬ 
stant  frequency,  and  that  the  term  cutoff  always  refers 
to  the  transverse  component  of  the  wave  vector,  k^,  not 
the  frequency,  o).  Isofrequency  contours,  a)(k)  =  const, 
show  the  required  relationship  between  k^  and  k^  for  plane 
wave  solutions  (Fig.  2). 

The  general  relationship  between  the  directions  of  en¬ 
ergy  and  phase  velocity  for  waves  propagating  within  an 
indefinite  medium  can  be  found  by  calculating  the  group 
velocity,  =  Vk6;(k).  Vg  specifies  the  direction  of  en¬ 
ergy  flow  for  the  plane  wave,  and  is  not  necessarily 
parallel  to  the  wave  vector.  Vkcu(k)  must  lie  normal  to 
the  isofrequency  contour,  m(k)  ==  const,  as  illustrated  in 
Fig.  2.  Calculation  of  Vkm(k)  from  the  dispersion  rela¬ 
tion,  Eq.  (3),  determines  which  of  the  two  possible  nor¬ 
mal  directions  yields  increasing  o)  and  is  thus  the  correct 
group  velocity  direction.  Performing  an  implicit  differ¬ 
entiation  of  Eq.  (3)  leads  to  a  result  for  the  gradient  that 
does  not  require  square  root  branch  selection,  removing 
any  sign  confusion.  To  obtain  physically  meaningful  re¬ 
sults,  a  causal,  dispersive  response  function,  must 
be  used  to  represent  the  negative  components  of  e  and 
/i,  since  these  components  are  necessarily  dispersive 
[14].  The  response  function  should  assume  the  desired 
(negative)  value  at  the  operating  frequency,  and  satisfy 
the  causality  requirement  that  d{^o))/do}  ^  1  [4,14]. 
Combining  this  with  the  derivative  of  Eq.  (3)  determines 
which  of  the  two  possible  normal  directions  applies, 
without  specifying  a  specific  functional  form  for  the 
response  function.  Figure  2  relates  the  direction  of  the 
group  velocity  to  a  given  material  property  tensor  sign 
structure. 

Having  calculated  the  energy  flow  direction,  we  can 
determine  the  refraction  behavior  of  indefinite  media  by 
applying  two  rules:  (i)  The  transverse  component  of  the 
wave  vector,  k^,  is  conserved  across  the  interface,  and 
(ii)  energy  carried  into  the  interface  from  free  space  must 
be  carried  away  from  the  interface  inside  the  media;  i.e., 
the  normal  component  of  the  group  velocity,  must 
have  the  same  sign  on  both  sides  of  the  interface.  Figure  2 
shows  typical  refraction  diagrams  for  the  three  types  of 
media  that  support  propagation. 

To  illustrate  the  unique  possibilities  associated  with 
indefinite  media,  we  recall  that  a  motivating  factor  in 
recent  metamaterials  research  has  been  the  propect  of 
near-field  focusing.  A  planar  slab  with  isotropic  s  ^  yu  = 
-1  can  act  as  a  lens  with  resolution  well  beyond  the 
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FIG.  2  (color).  Material  property  tensor  forms,  dispersion 
plot,  and  refraction  diagram  for  four  classes  of  media.  Each 
of  these  media  has  two  subtypes:  one  positive  (blue)  and  one 
negative  (red)  refracting,  with  the  exception  that  always  cutoff 
media  do  not  support  propagation  and  refraction.  The  disper¬ 
sion  plot  shows  the  relationship  between  the  components  of  the 
wave  vector  at  fixed  frequency,  (horizontal  axis)  is  always 
real,  (vertical  axis)  can  be  real  (solid  line)  or  imaginary 
(dashed  line).  (The  closed  contours  are  shown  circular,  but  can 
more  generally  be  elliptical.)  The  same  wave  vector  and  group 
velocity  vectors  are  shown  in  the  dispersion  plot  and  the 
refraction  diagram,  shows  direction  only.  The  shaded  diago¬ 
nal  tensor  elements  are  responsible  for  the  shown  behavior  for 
electric  y  polarization,  the  unshaded  diagonal  elements  for 
magnetic  y  polarization. 

diffraction  limit  [6].  It  is  difficult,  however,  to  realize 
significant  subwavelength  resolution  with  an  isotropic 
negative  index  material,  as  the  required  exponential 
growth  of  the  large  k^  field  components  across  the  nega¬ 


tive  index  lens  leads  to  extremely  large  field  ratios. 
Sensitivity  to  material  loss  and  other  factors  can  signifi¬ 
cantly  limit  the  subwavelength  resolution. 

By  combining  positive  and  negative  refracting  layers  of 
never  cutoff  indehnitQ  media,  we  can  produce  a  compen¬ 
sated  bilayer  that  accomplishes  near-field  focusing  in  a 
similar  manner  to  the  perfect  lens,  but  with  significant 
advantages.  Figure  2  indicates  that,  for  the  same  incident 
plane  wave,  the  z  component  of  the  transmitted  wave 
vector  is  of  opposite  sign  for  these  two  materials. 
Combining  appropriate  lengths  of  these  materials  results 
in  a  composite  indefinite  medium  with  unit  transfer  func¬ 
tion.  We  can  see  this  quantitatively  by  computing  the 
general  expression  for  the  transfer  function  of  a  bilayer 
using  standard  boundary  matching  techniques  [15], 

T  =  -  Zo)(l  +  Z,)(l  -  Zz) 

+  -  Zo)(l  -  Z,)(l  +  Zz) 

+  +  Zo)(l  -  Zi)(l  -  Z2) 

+  +  Zq)(i  +  z,)(l  +  Zz)]-‘.  (4) 


The  relative  effective  impedances  are  defined  as 


ry  —  ^zl 

^0 - r“» 


2  ^  fj'x\  ^z2 
Mjc2  ^zl 


where  k,  q^,  and  q2  are  the  wave  vectors  in  vacuum  and 
the  first  and  second  layers  of  the  bilayer,  respectively.  The 
individual  layer  phase  advance  angles  are  defined  as  ^  = 
qz\Lx  and  ij/  =  where  Lj  is  the  thickness  of  the  first 

layer  and  L2  is  the  thickness  of  the  second  layer.  If  the 
signs  of  and  q^2  opposite  as  mentioned  above,  the 
phase  advances  across  the  two  layers  can  be  made  equal 
and  opposite,  -h  =  0.  If  we  further  require  that  the 
two  layers  are  impedance  matched  to  each  other,  Zj  =  1, 
then  Eq.  (4)  reduces  to  7  =  1  (the  transfer  function  of 
free  space  is  7  =  ^*^^(^1+^2)),  in  the  absence  of  loss,  the 
material  properties  can  be  chosen  so  that  this  occurs  for 
all  values  of  the  transverse  wave  vector,  kj^.  Transfer 
functions  with  varying  amounts  of  loss  added  are  shown 
in  Fig.  3. 

A  proposed  implementation  is  shown  in  Fig.  1.  The 
elements  shown  in  the  top  and  bottom  of  the  figure  will 
implement  media  that  focuses  electric  y-polarized  and 
magnetic  y-polarized  waves,  respectively.  Combining 
the  two  structures  forms  a  bilayer  that  is  x~y  isotropic 
due  to  the  symmetry  of  the  combined  lattice.  This  sym¬ 
metry  and  the  property  jjl  =  e  yield  polarization  inde¬ 
pendence.  The  materials  are  formed  from  split  ring 
resonators  and  wires  with  numerically  and  experimen¬ 
tally  confirmed  effective  material  properties  [5].  Each 
split  ring  resonator  orientation  implements  negative  per¬ 
meability  along  a  single  axis,  as  does  each  wire  orienta¬ 
tion  for  negative  permittivity. 

While  compensated  bilayers  of  indefinite  media  ex¬ 
hibit  reduced  impedance  mismatch  to  free  space  and 
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FIG.  3.  Magnitude  of  the  transfer  function  vs  transverse  wave 
vector,  for  a  bilayer  composed  of  positive  and  negative 
refracting  never  cutojf  Material  property  elements  are  of 

unit  magnitude  and  layers  of  equal  thickness,  A.  A  loss  pro¬ 
ducing  imaginary  part  has  been  added  to  each  diagonal  com¬ 
ponent  of  E  and  jjL,  with  values  0.001,  0.002,  0.005,  0.01,  0.02, 
0.05,  and  0.1  for  the  darkest  to  the  lightest  curve.  For  compari¬ 
son,  a  single  layer,  isotropic  near  field  lens  is  shown  dashed 
The  single  layer  has  thickness  A,  and  e  =  /x,  =  -1  +  0.001/. 


high  transmission,  uncompensated  semi-infinite  indefi¬ 
nite  media  exhibit  unique  high  reflection  properties. 
Aside  from  cutojf  materials,  other  classes  of  indefinite 
media  have  a  reflection  coefficient  amplitude  near  unity 
for  incident  propagating  waves.  The  phase  of  the  reflec¬ 
tion  coefficient,  however,  varies,  as  illustrated  here  for 
positive  refracting  anti-cutojf  media.  The  reflection  co¬ 
efficient  for  electric  y  polarization  is  given  by 


f^x^z 


(6) 


where  k  and  q  are  the  wave  vectors  in  vacuum  and  the 
media,  respectively.  For  unit  magnitude  anti-cutojf  mei~ 
terial  we  have,  from  Eq,  (3), 

q\  = -^  + -k\.  (7) 

Thus,  The  correct  sign  for  positive  refracting 

media,  +,  is  determined  by  the  requirement  that  the  fields 
must  not  diverge  in  the  domain  of  the  solution.  We  then 
have 


The  magnitude  of  the  reflection  coefficient  is  unity,  with  a 
phase  of  '-90°  for  propagating  modes  of  all  incident 
angles.  An  electric  dipole  antenna  placed  A/8  away 
from  the  surface  would  thus  be  enhanced  by  interaction 
with  this  “mirror”  surface.  Customized  reflecting  surfa¬ 
ces  are  of  practical  interest,  as  they  enhance  the  efficiency 
of  nearby  antennas,  while  at  the  same  time  providing 
shielding  [16,17].  Furthermore,  an  interface  between 
unit  cutojf  and  anti-cutojf  media  has  no  solutions  that 
are  simultaneously  evanescent  on  both  sides,  i.e.,  surface 
modes  (Fig.  2).  In  many  communications  applications. 


the  energy  lost  to  the  excitation  of  surface  modes  is 
undesirable,  as  it  represents  loss  of  signal  power. 

In  conclusion,  we  have  begun  to  explore  the  properties 
of  media  with  indefinite  e  and  pi  tensors.  Consideration  of 
layered  structures  has  led  to  useful  and  interesting  reflec¬ 
tion  and  refraction  behavior,  including  a  new  mechanism 
for  subdiffraction  focusing.  We  note  that  neither  the 
analysis  nor  the  fabrication  of  these  media  is  compli¬ 
cated,  and  thus  anticipate  other  researchers  will  quickly 
assimilate  the  principles  and  design  structures  with 
unique  and  technologically  relevant  properties. 

We  thank  Claudio  Parazzoli  (Phantom  Works,  Boeing) 
for  motivating  this  work.  This  work  was  supported 
by  DARPA  (Contract  No.  972-01-2-0016),  ONR  (Con¬ 
tract  No.  N00014-00- 1-0632),  and  AFOSR  (Contract 
No.  F49620-01- 1-0440). 

Note  added — ^While  preparing  this  manuscript,  we 
learned  that  a  focal  spot  of  A/25  was  experimentally 
achieved  in  a  lumped  circuit  analog  of  the  bilayer  geom¬ 
etry  [18]. 
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A  planar  slab  of  material,  for  which  both  the  permittivity  and  permeability  have  the  values  of  -I, 
can  bring  not  only  the  propagating  fields  associated  with  a  source  to  a  focus,  but  can  also  refocus 
the  nonpropagating  near  fields,  thereby  achieving  resolution  beyond  the  diffraction  limit.  We  study 
the  sensitivity  of  this  subwavelength  focus  to  the  slab  material  properties  and  periodicity,  and  note 
the  connection  to  slab  surface  plasmon  modes.  We  conclude  that  significant  subwavelength 
resolution  is  achievable  with  a  single  negative  index  slab,  but  only  over  a  restrictive  range  of 
parameters.  ©  2003  American  Institute  of  Physics,  [DOI:  10.1063/1.1554779] 


A  wave  incident  on  the  interface  between  two  materials 
whose  indices  of  refraction  are  of  opposite  sign  will  undergo 
negative  refraction.  This  recently  appreciated  propagation 
phenomenon  has  been  predicted  to  lead  to  numerous  inter¬ 
esting  optical  phenomena,  including  a  superfocusing  effect. 
Veselago  showed  theoretically  that  a  planar  slab  with  index 
of  n  =  -“  1  could  focus  the  rays  from  a  nearby  electromag¬ 
netic  source  to  an  image  on  the  opposite  side  of  the  slab.^  In 
a  more  recent  analysis,^  it  was  shown  that  in  addition  to  the 
far-field  components  associated  with  a  source  being  brought 
to  a  focus  by  the  slab,  the  nonpropagating  near-field  compo¬ 
nents  could  also  be  recovered  in  the  image.  It  was,  therefore, 
proposed  that  the  image  created  by  a  planar  slab  could,  in 
principle,  contain  all  of  the  information  associated  with  the 
source  object,  thereby  achieving  resolution  well  beyond  that 
of  the  diffraction  limit.  For  this  reason,  the  slab  was  de¬ 
scribed  as  a  perfect  lens.  We  maintain  this  description  here, 
referring  specifically  to  a  planar  slab  of  continuous  material 
with  =  - 1  and  e  =  - 1  (no  losses)  as  a  perfect  lens. 

The  resolution  enhancement  associated  with  the  perfect 
lens  was  a  surprising  result,  stimulated  by  the  experimental 
demonstration^  of  a  left-handed  material  at  microwave  fre¬ 
quencies,  for  which  e<0  and  )a<0.  However,  far  from  be¬ 
ing  a  continuous  material,  the  measured  sample  was  com¬ 
prised  of  two  interlaced  periodic  arrays  of  copper  elements, 
one  array  being  composed  of  split  ring  resonators,^  and  the 
other,  wires.^ 

The  values  of  the  electromagnetic  parameters  and  the 
spatial  periodicity  render  the  experimental  sample  distinct 
from  the  idealized  perfect  lens.  The  question  then  arises  as  to 
whether  or  not  focusing  beyond  the  diffraction  limit  can  ac¬ 
tually  be  observed  using  any  practically  obtained  or  fabri¬ 
cated  material,  or  even  can  be  simulated  using  standard  nu¬ 
merical  methods  (e.g.,  finite  difference  or  finite  element) 
which  inevitably  approximate  the  ideal  situation.  We  explore 
the  inherent  limitations  associated  with  realizable  materials, 
and  the  expected  impact  on  the  focusing  properties  of  a  slab. 

The  planar  geometry  we  consider  here  allows  a  straight- 


“^Electronic  mail:  drs@sdss.ucsd.edu 


forward  analysis  to  be  implemented,  as  was  used  in  Ref.  2. 
The  fields  from  an  arbitrary  electromagnetic  source  are  ex¬ 
panded  in  a  Fourier  series  over  homogeneous  and  inhomo¬ 
geneous  plane  waves.  The  influence  of  the  slab  on  each  plane 
wave  component  can  easily  be  determined  by  a  standard 
transfer  matrix  technique.  Restricting  the  field  variation  to 
one  transverse  direction,  with  the  electric  field  having  S  po¬ 
larization,  we  find  the  field  expansion 

E(jc,z,r)  =  X  E{kx)txp{ik^Z‘^ikjpc-io}t)y,  (1) 

where  k^  and  kj^  are  the  components  of  the  wave  vector  nor¬ 
mal  to  and  parallel  to  the  slab,  respectively.  Outside  the  slab, 
the  wave  equation  leads  to  the  usual  dispersion  relation  re¬ 
lating  the  frequency  o)  and  the  components  of  the  wave  vec¬ 
tor,  or  k^-  ^J(o^/c^-kl.  Those  modes  for  which  \kx\<oy/c 
are  propagating,  while  those  for  which  \kx\>o)/c  decay  eva- 
nescently  along  the  propagation  direction  (z).  As  was  shown 
in  Ref.  2,  it  is  the  latter  inhomogeneous  modes  that  are  re¬ 
sponsible  for  image  resolution  beyond  the  diffraction  limit. 
The  ability  to  recover  these  components  in  an  image  is  what 
distinguishes  a  negative  index  structure  from  all  other  posi¬ 
tive  index  materials. 

For  each  plane  wave  component,  we  can  determine  a 
transfer  function,  defined  as  the  ratio  of  the  field  at  the  image 
plane  to  that  at  the  object  plane.  For  an  5-polarized  plane 
wave  incident  on  a  slab  of  thickness  d  and  arbitrary  values  of 
e  and  fJL,  the  transfer  function  has  the  form 


where  we  have  now  defined  q^=  yjk^— epof^Tc^,  and  rede¬ 
fined  kz=y[k^-^oPT^.  The  expression  for  P-polarized 
waves  is  similar  to  Eq.  (2),  with  the  explicitly  appearing  pi 
replaced  by  e.  We  apply  all  of  our  arguments  to  the 
iS-polarization  terms  in  this  letter,  as  the  results  for  P  polar¬ 
ization  follow  trivially. 
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FIG.  1-  Transfer  function  for  a  left-handed  slab.  For  the  perfect  lens,  the 
transfer  function  would  be  unity  for  all  /Cj .  However,  the  deviation  of  either 
the  real  or  imaginary  part  of  ji  limits  the  range  of  ,  so  that  the  slab  acts  as 
a  low-pass  filter.  The  losses,  inherent  to  left-handed  materials,  also  remove 
the  singularities  that  appear  in  the  transfer  function. 


When  fx  is  positive,  the  first  term  in  brackets  will  domi¬ 
nate  the  behavior  of  the  transmitted  wave  for  sufficiently 
large  d  and  \kj^\>(ofc,  with  the  evanescent  components  de¬ 
caying  exponentially  through  the  slab.  When  both  s  and  fx 
are  equal  to  —1,  however,  the  normally  dominant  solution 
vanishes,  and  I  for  every  component  from  the  source 
field — homogeneous  or  inhomogeneous,  thus,  exactly  repro¬ 
ducing  the  source  field  in  the  image  plane.  But  the  balance  is 
delicate;  any  deviation  from  the  perfect  lens  condition,  how¬ 
ever  small,  will  result  in  an  imperfect  image  that  degrades 
exponentially  with  slab  thickness  d,  until  the  usual  diffrac¬ 
tion  limit  is  reached.  This  sensitivity  has  been  noted  by  other 
researchers.^’^ 

The  effect  on  the  image  of  deviations  from  the  perfect 
lens  condition  for  a  given  slab  thickness  can  be  estimated  by 
determining  those  values  of  parameters  which  cause  the  two 
terms  in  Eq.  (2)  in  brackets  to  be  roughly  equivalent.  For 
small  deviations  from  the  perfect  lens  condition,  we  can  find 
an  approximate  expression  for  the  resolution  of  the  slab.  For 
e  =  —  1  and  /i=  —  1  +  Sfx,  the  two  terms  in  the  denominator 
of  Eq.  (2)  are  of  approximately  the  same  magnitude  when 
2k^d=  -  \xi\{8ixl2)^\.  Assuming  that  this  limit  occurs  when 
the  value  of  is  large,  we  have  we  can 

replace  with  .  At  the  maximum  k^ ,  the  minimum  re¬ 
solvable  feature  will  be  on  the  scale  of 
thus,  find  the  resolution  enhancement,  /?  =  A:jp/A:o»  of  the  lens 
as  a  function  of  small  deviations  in  the  permeability  is 
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(3) 


For  example,  if  1.5,  we  find  that  to  achieve  an  R  of  10, 
Sfx  must  be  no  greater  than  ~6X  10”^^!  However,  for  X/J 
=  10,  Sfx  can  vary  by  as  much  as  ~0.002  to  achieve  the  same 
resolution  enhancement.  Deviations  in  either  the  real  or  the 
imaginaiy  part  of  the  permeability  will  result  in  the  same 
resolution.  The  effect  of  varying  e  (for  S  polarization)  is 
much  smaller  than  that  associated  with  /x,  and  we  do  not 
discuss  the  effects  of  this  variation  further. 

The  transfer  function  has  poles  that  occur  when: 


qz 


or 


qz  ‘ 


(4) 


Equation  (4)  corresponds  to  the  dispersion  relations  for  slab 
plasmon  polaritons.^’^  For  certain  values  of  kj^ ,  these  disper¬ 
sion  relations  are  satisfied,  and  resonant  surface  modes  can 
exist  on  the  slab.  The  direct  excitation  of  these  surface 
modes  for  imaging  applications  is  undesirable,  as  the  corre¬ 
sponding  kx  components  will  be  disproportionately  repre¬ 
sented  in  the  image.  Yet,  the  existence  of  these  resonances  is 
essential,  as  the  recovery  of  the  evanescent  modes  can  be 
seen  as  the  result  of  driving  the  surface  plasmon  far  off  reso¬ 
nance.  The  appearance  of  resonances  can  be  seen  in  Fig.  (1), 
which  shows  the  transfer  function  as  a  function  of  k^ . 

While  losses  severely  limit  the  obtainable  resolution  of 
the  focus  for  a  negative  index  slab,  periodicity,  as  exists  for 
example  in  structured  metamaterials^  or  in  the  sampling  as¬ 
sociated  with  finite-difference  simulations,  also  imposes  a 
significant  resolution  limitation.  The  effect  of  periodicity  is 
minimal  for  propagating  plane  waves  whose  wavelength  is 
much  larger  than  the  repeated  unit-cell  size;  however,  peri¬ 
odicity  has  a  significant  effect  on  the  recovery  of  the  non¬ 
propagating  components  having  large  transverse  wave  num¬ 
ber. 

We  can  introduce  periodicity  into  this  analysis  in  an  ap¬ 
proximate  manner.  The  wave  equation  in  a  medium,  assum¬ 
ing  that  E  is  polarized  in  the  y  direction,  is 


We  have  assumed  that  the  periodic  variation  in  the  index  is 
only  in  the  transverse  (a:)  direction.  Under  these  assumptions, 
n^(x)  =  n^(x-\-a)j  where  a  is  the  repeat  distance. 

Due  to  the  periodicity  in  n^ix),  Eq.  (5)  in  Fourier  space 
will  involve  sums  over  reciprocal  lattice  vectors 
=  2n7rfa.  However,  to  obtain  a  rough  estimate  of  the  limi¬ 
tation  that  periodicity  imposes  on  the  resolution  enhance¬ 
ment,  we  solve  Eq.  (5)  for  an  index  having  the  form 


n^(jc)=  1  +  2  A  sin 


27rx 

a 


(6) 


The  validity  of  this  approximate  expression  can  be  seen  by 
comparison  with  Fig.  1.  For  example,  for  a  deviation  of 
0.005  and  \/d=  10,  the  numerically  computed  transfer 
function  in  Fig.  1  shows  that  the  range  of  k^  values  near 
unity  is  approximately  up  to  Eq*  (3)  predicts  R 

-9.5, 

The  dependence  of  the  resolution  enhancement  R  on  the 
deviation  from  the  perfect  lens  condition  is  critical.  The  ratio 
k/d  (wavelength  to  slab  length)  dominates  the  resolution,  the 

logarithm  term  being  a  relatively  weakly  varying  function. 
Downloaded  18  Apr  2003  to  132.239.27.179.  Redistribution  subji 


Using  this  form  in  Eq.  (5)  and  assuming  the  periodic  modu¬ 
lation  is  sufficiently  weak  that  only  two  bands  need  be  con¬ 
sidered,  we  find  the  modified  dispersion  relation  by  evaluat¬ 
ing  the  determinant 

^0  >^x  ‘Iz  ^  ,  , 

where  g  =  2'Trla  and  kQ  —  oylc.  When  A<^1,  the  bands  only 

weakly  couple,  and  we  ignore  the  effects  of  higher  mode 
to  AlP  license  or  copyright,  see  http://o)ps.aip,org/aplo/aplcr.jsp 
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excitation  to  find  an  expression  for  ,  When  the  two  terms 
in  the  denominator  of  Eq.  (2)  are  roughly  equal,  we  then 
have  the  condition 


9,d=-ln 


1 


2  ki 


+  sinh 


-1 


8^-K 


ikW 


(8) 


In  the  limit  kQLIg<\,  the  resolution  enhancement  is  ap¬ 
proximately 


X  1  /  xM  X 

d- 


(9) 


A  similar  resolution  enhancement  limit  to  that  obtained 
by  a  variation  in  material  parameter  is  thus  imposed  by  the 
introduction  of  periodic  modulation.  Similar  to  Eq.  (6),  the 
effects  of  periodicity  enter  logarithmically,  and  again  lead  to 
a  critical  dependence.  As  an  example,  an  of  10  for  a  slab 
with  X/^=10  and  A^-'l  requires  a  periodicity  of  ~X/20. 
While  admittedly  of  limited  quantitative  use,  Eq.  (9)  shows 
that  the  inherent  periodicity  in  metamaterials  will  impose  a 
limitation  on  the  resolution  of  the  lens.  This  same  limitation 
will  result  also  from  numerical  methods  that  model  a  con¬ 
tinuous  material  by  evaluating  the  fields  at  a  finite  number  of 
sampling  points  periodically  spaced.  This  limitation  has  un¬ 
doubtedly  complicated  numerical  attempts  to  observe  the  su- 
perfocusing  effect.*®  A  more  recent  numerical  study  has  con¬ 
cluded  that,  within  the  parameter  range  determined  from  Eq, 
(6)  and  with  a  fine  enough  discretization  grid,  enhanced  reso¬ 
lution  can  indeed  be  observed  in  finite-difference 
calculations.** 

It  should  be  noted  that  we  have  applied  a  definition  for 
resolution  in  this  work  that  may  not  be  appropriate  for  all 


applications.  For  instance,  calculations  (not  shown)  suggest 
that  the  excitation  of  surface  slab  plasmons  leads  to  a  broad 
background  in  the  image  plane,  with  the  subwavelength  fea¬ 
tures  superposed.  However,  for  certain  applications,  the  sub¬ 
wavelength  information  obtained  from  such  an  image  may 
be  of  value. 

Our  conclusion  from  this  analysis  indicates  that  the  per¬ 
fect  lens  effect  exists  for  a  fairly  restricted  region  of  param¬ 
eter  space.  Yet,  demanding  as  these  specifications  are, 
achieving  subwavelength  resolution  is  possible  with  current 
technologies.  Negative  refractive  indices  have  been  demon¬ 
strated  in  structured  metamaterials,  and  such  materials  can  be 
engineered  to  have  tunable  material  parameters  so  as  to 
achieve  the  optimal  conditions.  Losses  can  be  minimized  in 
structures  utilizing  superconducting  or  active  elements. 
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